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ABSTRACT 

The observatio n of a broad excess of sub-TeV cosmic rays compatible with the direction of the 
heliospheric tail (jNagashima et al]|1998f) and the discovery o f two significant lo calized excess regions 
of multi-TeV cosmic rays by the MILAGRO collaboration ([Abdo et al.l 120081 ). also from the same 
region of the sky, have raised questions on their origin. In particular, the coincidence of the most 
significant localized region with the direction of the heliospheric tail and the small angular scale of the 
observed anisotropy (~ 10°) is suggestive a local origin and of a possible connection to the low energy 
broad excess. Cosmic ray acceleration from magnetic reconnection in the magnetotail is proposed as 
a possible source of the energetic particles. 

Subject headings: magnetic fields- MHD- solar wind-acceleration of particles-cosmic rays 



1. INTRODUCTION 

It is known that cosmic rays arrival direction has an 
energy dependent large angular scale anisotropy with an 
amplitude of order 10~ 4 — 10~ 3 . The first comprehen- 
sive observation of this anisotropy was provided by a 
network of muon telescopes sensitive to sub-TeV ener- 
gies a nd located at different latitudes (iNagashima et all 
1998). More recently, an anisotropy was also ob- 
served in the multi-TeV energy range by the Tibet 
AS7 array (lAmenomori et alJl200a l. Super- Kamiokande 
(Guill ian et al.l 120071 ) and bv MILAGRO (|Abdo et al.l 
2009), and the first high statistics observation in the 
southern hemisphere in the 10 TeV reg ion, is being re- 
ported by IceCube (|Abbasi et all 120101) . The origin of 
the large angular scale anisotropy in the cosmic rays ar- 
rival direction is still unknown. The structure of the local 
interstellar magnetic field is likely to have an important 
role. However the combined study of the anisotropy en- 
ergy and angular dependency, its time modulation and 
angular scale structure seem to suggest that the obser- 
vation might be a combination of multiple superimposed 
effects, caused by phenomenologies at different distances 
from Earth. 

In this context, particular interest is derived from the 
observation of a broad excess of sub- TeV cosmic rays 
in a portion of the sky compatible with the direction 
of th e heliospheric ta il (or heliotail) (Naga shlma et all 
119981: lHall et allfl999t ) (see g2]). The heliotail is the re- 
gion of the heliosphere downstream the interstellar mat- 
ter wind delimited within the heliopause, i.e. the bound- 
ary that separates the solar wind and interstellar plasmas 
(|Izmodenov fc Kallenbach 2006). The observed excess 
was attributed to some unknown anisotropic process con- 
nected with the heliotail (thus called tail- in excess). The 
gyro-radius of sub- TeV cosmic protons is less than about 
200 AU (in a ~1 fiG interstellar magnetic field), which 
is approximately the size of the heliosphere and, most 
likely, smaller than the width and length of the heliotail. 



The persistence of the cosmic ray anisotropy structure in 
the multi-TeV energy range makes it challenging to link 
this observation to the heliosphere. Although the un- 
known size and extension of the heliotail contributes to 
the uncertainty on the energy scale at which heliospheric 
influence on cosmic rays starts to be negligible. However, 
we know that the observations of multi-TeV cosmic rays 
anisotropy show small angular scale patterns superim- 
posed to the smooth broad structure of the tail-in excess, 
which is suggestive of a local origin, i.e. within the helio- 
tail. With the same technique used in gamma ray detec- 
tion to estimate the background and search for sources 
of gamma rays, the MILAGRO collaboration discovered 
two localized excess regions in the cosm ic rays arrival di- 
rection distribution ()Abdo et al.ll2008() . The same excess 
regi ons were reported by the ARGO-YBJ air shower ar- 
ray (|Vernetto et al.ll2009l ). The strongest and more local- 
ized of them (with an angular size of about 10°) coincides 
with the direction of the heliotail. The peculiarity of 
such an observation triggered an astrophysical interpre- 
tation based on the possibility that cosmic rays acceler- 
ated by the supernova that produced Geminga pulsar are 
focus sed by an ad-hoc inter s tellar magnetic field struc- 
ture (jSalvati fc Saccol 120081 : iDrurv fc Aharonianl 120081 ) 
(see ©. 

The localized regions lie in the same portion of the sky 
that is dominated by the broad tail-in excess at lower 
energy, and, although it might be coincidental, we inter- 
pret this as manifestations of the same phenomenology 
at different energies. 

It is proposed that both sub- TeV tail-in excess and 
the multi-TeV localized excess of cosmic rays might be 
caused by magnetic reconnection in the heliosphere and, 
in particular in the heliotail, where the distance scale 
might be long enough to induce sufficient acceleration 
at high energy. The very idea of appealing to magnetic 
reconnection for the acceleration of energetic particles 
can be traced back to pioneeding works by Giovanelli 
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(1946) and Dungey (1953). The uncertainties with un- 
derstanding of fast reconnection were one of the imped- 
iments for applying the process to energetic particle ac- 
celeration (see Lazarian & Opher 2009). We appeal to 
the model of reconnection of weakly stochastic field in 
Lazarian & Vishniac (1999), which was identified as a 
cause of First Order Fermi acceleration (see de Gouveia 
dal Pino & Lazarian 2005, Lazarian 2005). 

In what follows we present the observational evidence 
for the existence of the cosmic ray excess in the direction 
of the solar system magnetotail in fj2j discuss existing 
explanations of this excess in Ej3j The structure of the 
magnatotail with magnetic field reversals arising from 
the solar cycle is presented in §4] and the mechanism of 
acceleration of cosmic rays in the magnetotail is outlined 
in The discussion of the results and a short summary 
are given by [J6] and $3 respectively. 

2. THE OBSERVATIONS 

The observation of t he large angular scale anisotropy of 
sub-T eV cosmic rays (|Nagashima et alJll998t lHall et al.l 
1999) revealed the evidence of a superposition of two 
different modulations in arrival direction. One with a 
sidereal variation identified with an extended deficit cen- 
tered around 12 hours that seems to extend mostly across 
the northern hemisphere (the so-called loss cone). And 
one with a sidereal variation identified with a broad ex- 
cess centered around 6 hours, with half opening angle of 
about 68° that comprises the direction of the heliotail, 
and extended across part of the northern and the south- 
ern hemispheres (the tail-in excess) . Figure [1] shows the 
combined observations of the anisotropy of sub-TeV cos- 
mic rays from telescopes at different latitudes. 

The global anisotropy amplitude is found to increase 
with energy up to about 5-10 TeV, however while the 
loss-cone structure seems to maintain a similar shape up 
to the multi-TeV range, the tail-in excess is still some- 
what persistent in the multi-TeV range, but its broad 
struct ure appears to dissolve to smaller angular scale 
spots (|Amenomori et abl 120061 ). The apparent seasonal 
modulation of the tail-in excess, with a minimum ampli- 
tude in summer and a maximum (a factor of four larger) 
in winter, provides a compelling connection to the helio- 
tail. 

Figure [2] shows the multi-TeV cosmic ray arrival direc- 
tion map, from the MILAGRO collaboration, obtained 
by eliminating anisotropics with angular structures wider 
than ^30°. The small scale structure is evidenced in this 
map and it shows two highly significant (more than 12 
cr) localized excess regions in the cosmic rays arrival di- 
rection. Both regions are inconsistent with gamma ray 
emission with high confidence and therefore are claimed 
to be dominated by cosmic rays. They are found to have 
a constant yearly excess over the seven year period of col- 
lected data, however both of them were lowest in summer 
and highest in winter, with a \ 2 probability relative to 
a constant fractional excess of only 5% in each region. 
The strongest and more localized of them (called region 
A, with a fractional excess of ~ 6 x 10~ 4 ) coincides with 
the direction of the heliotail (the black dot in Figure [21 
with right ascension a ps 74° and declination S *=s +17° in 
equatorial coordinates) . The corresponding energy spec- 
trum was de-convoluted using the energy-dependent ex- 
perimental observables. Figure [3] shows the result of the 
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Fig. 1. — Cosmic ray sidereal daily variation S(t) from differ- 
ent muon telescope station or air shower detectors : from top 
to bottom, Mount Norikura air shower array (F3, shown for 
reference at high energy and with the amplitude multiplied 
by 1/4 to account for the larger anisotropy amplitude at high 
energy), Nagoya station (looking 30° north- ward, and look- 
ing vertically up- ward), Sakashita station (looking vertically 
up-ward and 41° south-ward) and Hobart station (looking 
vertically up-ward). Each station reports the corresponding 
latitude of the directional detection and the cosmic ray me- 
dian energy it is sensitive to. In order to clearly show the 
peaks and valleys, the 24 hour variation is repeatedly shown 
in a 2 days time interval. The error expresses the dispersion 



of the hourly relative intensity. FromlNagashima et al 
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Fig. 2. — Map of statistical significances from MILAGRO. 
A 10° bin was used to smooth the data, and the color scale 
gives the significance. The solid line indicates the galactic 
plane and the dash lines the galactic latitudes. The black dot 
indicates the heliotail. The fractional excess of region A is 
~_6x 10~ 4 , and of region B is ~ 4 x 10~ 4 . From A bdo et al.l 
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X 2 fit to the excess in region A assuming a pure proton 

_ E 

spectrum of the form E 1 ■ e E <= , where 7 is the cosmic 
ray spectral index and E c the cut-off energy. 
Figure 3 testifies that the estimated cosmic ray energy 
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Fig. 3. — Result o f a y 2 fit to the ex cess in region A (from 
top of Figure 4 in (|Abdo et all 12001 0. The la, 2a and 3cr 
allowed regions of the spectral index 7 and cut-off energy 
E c are indicated by the shaded regions (from black to gray, 
respectively) . 

spectrum is consistent to harder than the isotropic flux 
(at a 4.6 a level) with a cutoff, and with the most sign if- 
icant excess in the multi-TeV range (lAbdo et al.l l2008). 

The localized excess regions of multi-TeV cosmic rays 
cover the same portion of the sky where the tail-in excess 
was observed at cosmic ray energies below TeV, and have 
similar seasonal modulations. It is likely that these are 
two manifestations of the same phenomenology and that 
the heliotail has an important role. 

IceCube has recently reported a first view of the multi- 
TeV medium scale anisotropy in the southern hemi- 
sphere (i.e. only mo dulations that are smaller than 60°) 
(jDesiati et al.l l2010), which might add novel information 
to this observation. 

3. OTHER INTERPRETATIONS 

While no explanation has ever being attempted to ex- 
plain the broad sub- TeV tail-in excess, a number of inter- 
pretations have been provided to address the existence of 
the most significant localized excess of multi-TeV cosmic 
rays. 

Some proposed models rely on astrophysical origin. 
In ISalvati fc Saccol (|2008D it is noted that the two lo- 
calized excess regions observed by MILAGRO surround 
the present day apparent location of Geminga pulsar. 
The supernova that gave birth to the pulsar exploded 
about 340,000 years ago, when its distance to the Sun 
was estimated to be about 90 pc. Even if the proper 
motion of the pulsar induced by the explosion moved it 
further away (the present distance of Geminga pulsar is 
estimated to be about 155±35 pc), 10 TeV cosmic rays 
produced by the supernova have propagated about 65 
pc away, if we assume Bohm diffusiorQ from the source 
to here : approximately consistent with the distance of 
Geminga at the time of explosion. From this distance 
a total cosmic ray energy of about 1.5 ■ 10 49 erg must 
be emitted by the supernova to produce the observed 
fractional excess. This value is consistent with the com- 
monly required efficiency (~ 1%) with which a supernova 
energy output must be converted into cosmic rays if they 
are to maintain the galactic cosmic ray density. 

1 Bohm diffusion assumes scattering of a cosmic ray for every 
gyration (see Parker 1979). 



The major problem with this explanation is that Bohm 
diffusion through such large distances cannot possibly ex- 
plain the localized nature of the observed excesses, but 
it would rather produce at most a broad faint dipolar 
anisotropy in arrival direction. In addition, if Bohm scal- 
ing for the scattering in the immediate vicinity of the 
supernova shock is plausible, it seems unlikely that this 
regime persists during the propaga tion of cosmic rays 
throu gh the interstellar medium (|Drurv fc Aharonianl 
2008). On the other hand the structure of the inter- 
stellar magnetic field would very hardly maintain multi- 
TeV cosmic rays focussed within a 10° beam. The sup- 
posed opposite scenario of a free-stre aming of cosmic rays 
along a sort of magnetic "nozzle" (jDrurv fc Aharonianl 
2008), that would explain the localized nature of the ob- 
servation, would also be extremely unlikely also because 
the propagation would have been so fast that we would 
not have the observation in the first place anymore. It 
is possible to argue a scenario with a combination of 
slow Bohm diffusion regime (close to the supernova due 
to turbulence induced by the explosion) and fast free- 
streaming along magnetic field through the interstellar 
medium (which could partially explain the localized na- 
tures of the observed excesses). But this interpretation 
would have problems, as the time of propagation from the 
source should be made long, which contradicts the idea 
of localization of the intensive scattering only near the 
source. Perhaps some sort of leaky magnetic field bottles 
are formed near the source, which make the propagation 
slow compared to the Bohm diffusion time, thus mitigat- 
ing the problem. However this possibility would require 
fine tuning and we have not seen this idea discussed in 
the literature. 

The coincidence of the most significant localized ex- 
cess observed by MILAGRO with the heliotail, supports 
the idea that the heliosphere could somehow have a role. 
The possibility that we are seeing the effects of neutron 
production in the gravitationally fo cussed tail of the in- 
terste llar material was considered bv lDrurv fc A haronian 
(2008). As the Solar system surrounded by Solar Wind 
moves through the interstellar medium, the complex in- 
teraction between the two media create the heliotail. 
Cosmic rays propagating through the direction of the 
tail interact with the matter and magnetic fields to pro- 
duce neutrons and hence a localized excess of cosmic ray 
in that direction. But while the target size has about 
the right size compared to the decay length of multi-TeV 
neutrons (~ 0.1 pc), the increase of the gravitating mat- 
ter density is too low to account for the observed excess. 

While it is possible to argue that the large angular 
scale anisotropy of cosmic rays arrival direction might 
be generated by a combination of astrophysical phenom- 
ena, such a s the distribution of nearby recent supernova 
explosions (jErlvkin fc Wolfendald[2Tj06f) . propagation ef- 
fects (Battaner, Castellano & Masip 2009, Malkov et al. 
2010) and the structure of the interstellar magnetic field, 
it is more likely that small angular scale anisotropics are 
generated by some local effect. 

4. MAGNETIC FIELD STRUCTURE AT THE HELIO-TAIL 

Figure |4] represents the possible structure of the he- 
liotail which arises from the solar magnetic field cycles 
(Parker 1979). The magnetic fields of the opposite polar- 
ities emerge as the result of 11 year solar dynamo cycle. 
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Fig. 4. — A meridional view of the boundary sectors of the 
heliospheric currenty sheet and how the opposite sectors get tighter 
closer to the heliopause and into the heliotail. The thickness of the 
outflow regions in the reconnection region depends on the level of 
turbulence. The length of the outflow regions L depends on the 
mean geometry of magnetic field and turbulence Adapted from 
Nerney et al. 1995 and also Lazarian & Opher 2009. 

As the magnetic field is carried away by solar wind, the 
reversed magnetic field regions get accumulated in the 
magnetotail region. This is where reconnection is ex- 
pected to occur. 

Naturally, the actual heliotail is going to be turbulent, 
which is not represented by the idealized drawing in Fig- 
ure |H As the Alfven speed is smaller than the Solar wind 
speed, magnetic reconnection does not change the overall 
magnetic field structure. Nevertheless, as we discuss in 
§5, the effects of turbulence are very important from the 
point of view of magnetic reconnection and the particle 
acceleration that it entails. 

The simulations of the magnetotail are extremely chal- 
lenging (see Pogorelov et al. 2009ab) and have not been 
done with the sufficient resolution and extent. While we 
believe that future research will provide details necessary 
for quantitative modeling, the schematic representation 
of the magnetotail structure depicted in Figure 2] is true 
in terms of major features. In what follows, it will be 
used for describing the scenario for the origin of the cos- 
mic ray excess that we advocate in this paper. 

5. MAGNETIC RECONNECTION AND COSMIC RAY 
ACCELERATION 

Astrophysical plasmas are often highly ionized and 
highly magnetized (Parker 1970). The evolution of the 
magnetic field in a highly conducting fluid can be de- 
scribed by a simple version of the induction equation 



dB 
~dt 



V x (v x B - r]\7 x B 



(1) 



where B is the magnetic field, v is the velocity field, and 
rj is the resistivity coefficient. Under most circumstances 
this is adequate for discussing the evolution of magnetic 
field in an astrophysical plasma. When the dissipative 
term on the right hand side is small, as is implied by 
simple dimensional estimates, the magnetic flux through 
any fluid element is constant in time and the field topol- 
ogy is an invariant of motion. On the other hand, recon- 
nection is observed in the solar corona and chromosphere 
(Innes et al. 1997, Yokoyama & Shibata 1995, Masuda 
et al. 1994, Ciaravella & Raymond 2008), its presence 
is required to explain dynamo action in stars and galac- 
tic disks (Parker 1970, 1993), and the violent relaxation 
of magnetic fields following a change in topology is a 
promising process for the First order Fermi acceleration 




Sweet-Parker model 





Turbulent model 



blow up 



Fig. 5. — Upper plot: Sweet-Parker model of reconnection. The 
outflow is limited by a thin slot A, which is determined by Ohmic 
diffusivity. The other scale is an astrophysical scale L> A. Middle 
plot: Reconnection of weakly stochastic magnetic field according to 
LV99. The model that accounts for the stochasticity of magnetic 
field lines. The outflow is limited by the diffusion of magnetic 
field lines, which depends on field line stochasticity. Low plot: An 
individual small scale reconnection region. The reconnection over 
small patches of magnetic field determines the local reconnection 
rate. The global reconnection rate is substantially larger as many 
independent patches come together. The bottleneck for the process 
is given by magnetic field wandering and it gets comparable to L 
as the turbulence injection velocity approaches the Alfvenic one. 
From Lazarian, Vishniac & Cho (2004). 



of high energy particles in the universe (de Gouveia Dal 
Pino & Lazarian 2003, henceforth GL03, 2005, Lazarian 
2005, Drake et al. 2006, Lazarian & Opher 2009, Drake 
et al. 2010). Quantitative general estimates for the speed 
of reconnection start with two adjacent volumes with dif- 
ferent large scale magnetic field directions (Sweet 1958, 
Parker 1957). 

The speed of reconnection, i.e. the speed at which 
inflowing magnetic field is annihilated by Ohmic dissipa- 
tion, is roughly rj/ A, where A is the width of the tran- 
sition zone (see Figure [5]) . Since the entrained plasma 
follows the local field lines, and exits through the edges 
of the current sheet at roughly the Alfven speed, Va, the 
resulting reconnection speed is reduced compared to the 
Alfven speed by a large factor. This factor in the text- 
book Sweet-Parker model of reconnection is S 1 ^ 2 , where 
S = (LVa/t]) is the Lunquist number, where L is the 
length of the current sheet (Sweet 1958, Parker 1957, see 
also Parker 1979). 

In general, satisfying the conservation of mass condi- 
tion dictates that V rec ~ Va(L/A). Observations re- 
quire a speed close to Va, so this expression implies that 
L ~ A, i.e. that the region of over which magnetic flux 
tubes intersect should be comparable with the outflow re- 
gion. This can be achieved either via making L as small 
as the Ohmic diffusion region i.e. that the magnetic 
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field lines reconnect in an "X point" (Petscheck 1964), 
or the outflow region should get increased dramatically 
beyond the size that is predicted in the Sweet-Parker 
model. While for years the problem of fast reconnection 
was viewed as connected with proving of the stability of 
X point, the situation has changed recently with the sec- 
ond way of dramatically increasing the thickness of the 
outflow is becoming more popular and getting observa- 
tional support (see Ciaravella & Raymond 2008). 

The first model of X point reconnection was proposed 
by Petschek (1964). In this case the reconnection speed 
may have little or no dependence on the resistivitj0. The 
X point configuration is known to be unstable to collapse 
into a sheet in the MHD regime (see Biskamp 1996), but 
in a collisionless plasma it can be maintained through 
coupling to a dispersive plasma mode (Sturrock 1966). 
Recent years, have been marked by the progress in un- 
derstanding some of the key processes of reconnection 
in astrophysical plasmas. In particular, a substantial 
progress has been obtained by considering reconnection 
in the presence of Hall effect, which is described by the 
J x B term in Ohm's law: 

„ v x B J x B AirnJ , , 

E+ = —L- (2) 

c en e c cr 

where e is electron charge and n e is concentration of elec- 
trons. Numerical experiments showed that Hall-MHD 
reconnection is capable of supporting X-points and thus 
can make the reconnection fast, i.e. comparable to the 
Alfven speed (Shay et al. 1998, 2004). 

The condition at which Hall-MHD term gets impor- 
tant for the reconnection is that the ion skin depth 8i on 
is comparable with the Sweet-Parker diffusion scale A. 
The ion skin depth is a microscopic characteristic and it 
can be viewed at the gyroradius of an ion moving at the 
Alfven speed, i.e. Si on = Va/^ci, where w c ; is the cy- 
clotron frequency of an ion. In the heliotail for a proton 
we find that Si on ~ 10 3 km. Thus one can get the con- 
straint on the scale L for which Hall-MHD effects should 
dominate the reconnection: 

^«0. 2 (-^-) 1/2 ^ 4 <l, (3) 

where X m f p is the electron mean free path, where j3 p i is 
the ratio of thermal pressure to magnetic pressure (see 
more discussion in Yamada et al. 2006). We agrue in 
Appendix A that in realistic situations in turbulent me- 
dia the scales A || = L tur b (see Figure [5]) over which the 
microscale Sweet-Parker reconnection of individual tur- 
bulent patches may take place are much smaller than the 
scale of the system and therefore the collisionless effects 
take place within the heliotail. This, as we argue below, 
does not change the overall rates of magnetic reconnec- 
tion. 

A shortcoming of many discussions of magnetic re- 
connection is that the traditional setup does not in- 
clude ubiquitous pre-existing astrophysical turbulenc43 

2 In general, the reconnection is termed fast when the reconnec- 
tion velocity does not depend on the Lundquist number S or if it 
depends on ln(S'). In all other cases the large values of S make 
reconnection too slow for most of astrophysical applications. 

3 The set ups where instabilities play important role include 
Simizu et al. (2009a, b). For sufficiently large resolution of simula- 



(see Armstrong, Rickett & Spangler 1994, Elmegreen 
& Scalo 2004, McKee & Ostriker 2007, Lazarian 2009, 
Chepurnov & Lazarian 2010). As turbulence radically 
changes many astrophysical processes, the influence of 
turbulence on reconnection has attracted the attention 
of researchers for a long time (see Speizer 1970, Straus 
1988) . An extended discussion of turbulence role in accel- 
erating reconnection can be found in Mathaeus & Lamkin 
(1985, 1986). However, there the X-points created by 
turbulence together with the effects of compressibility 
and heating were identified as the means of accelerating 
reconnection. 

A very different approach to the effects of turbulence 
was adopted in Lazarian & Vishniac (1999, henceforth 
LV99). Their model does not appeal to any of the effects 
of turbulence-created X-points, compressibility or heat- 
ing and it is applicable to a wide range of astrophysical 
conditions. Fortunately, this approach provides a robust 
way of accelerating reconnection. Indeed, as we men- 
tioned above, the approach in LV99 is to consider ways 
to decouple the width of the plasma outflow region from 
the scale determined by Ohmic effects^ The plasma is 
constrained to move along magnetic field lines, but not 
necessarily in the direction of the mean magnetic field. 
In a turbulent medium the two are decoupled, and fluid 
elements that have some small initial separation will be 
separated by a large eddy scale or more after moving the 
length of the current sheet. As long as this separation is 
larger than the width of the current sheet, the result will 
not depend on 77. The mental picture presented in LV99 
is that the fluid follows magnetic field lines, which are 
not straight, but wander (see Figure [5] and also Lazarian 
et al. 2004, where this wandering was calculated numer- 
ically and compared with the analytical predictions in 
LV99). As a result, the thickness A of the fluid outflow 
is determined not by microphysical Ohmic diffusivity, but 
magnetic field wandering which for the injection velocity 
Vi of the order of Va is of the order of the turbulence 
injection scale I, i.e. Aw/. If the length of the current 
sheet L is of the order of I it is clear from Figure [5] that 
V rec can be comparable with Va- Note that LV99 con- 
siders generic 3D configurations of magnetic fluxes with 
non-zero magnetic guide, i.e. shared, magnetic field. The 
shared magnetic field is being ejected from the reconnec- 
tion region together with plasmas. The cosmic rays that 
we consider in this paper stay entrained on magnetic field 
lines. 

Two effects are the most important for understanding 
of the nature of reconnection in LV99. First of all, in 
three dimensions bundles of magnetic field lines can en- 
ter the reconnection region and reconnect there indepen- 
dently (see Figure [SJ , which is in contrast to two dimen- 
sional picture where in Sweet-Parker reconnection the 
process is artificially constrained. Then, the nature of 
magnetic field stochasticity and therefore magnetic field 
wandering (which determines the outflow thickness, as il- 

tions those set-ups are expected to demonstrate turbulence. Tur- 
bulence initiation is also expected in the presence of plasmoid ejec- 
tion (Shibata & Tanuma 200f ). Numerical viscosity constrains our 
ability to sustain turbulence via reconnection, however. 

4 In the Sweet-Parker reconnection both widths coincide and 
equal to A. In the LV99 model the outflow A is much wider than 
the thickness of the individual current sheets, while the Sweet- 
Parker is obtained as a degenerate case of no turbulence. 
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lustrated in Figure [5]) is very different in 2D and the real 
3D world (LV99). In other words, by removing artificial 
constraints on the dimensionality of the reconnection re- 
gion and the magnetic field being absolutely straight, 
LV99 explore the real-world astrophysical reconnection, 
which has become a topic of modern numerical studietO 
(see Dorelli & Bhattacharjee 2008, Kowal et al. 2009, 
Daughton et al. 2009). 

Analytical calculations in LV99 showed that the re- 
sulting reconnection rate is limited only by the width of 
the outflow region. That model predicts reconnection 
speeds close to the turbulent velocity in the fluid. More 
precisely, assuming isotropically driven turbulence char- 
acterized by an injection scale, I, smaller than the current 
sheet length L, LV99 obtained 

V rec « V A (l/L) 1/2 (Vi/Va) 2 , (4) 

where the turbulent injection velocity Vi is assumed to be 
less than Va- If L < I, the first factor in Eq. (j4} should 
be changed to {L/l) 1 / 2 (LV99). Taking into account that 
Vturb = Va(Vi/Va) 2 is the velocity at the scale where 
turbulence transits from weak MHD turbulence to the 
strong MHD turbulenc^B (see LV99, Lazarian 2006), one 
can rewrite Eq. (|4]) in the following way: 

Vrec « Vturb {l/L) 1/2 , (5) 

Note, that here "strong" means only that the eddies de- 
cay through nonlinear interactions in an eddy turn over 
time (see more discussion of the LV99). All the mo- 
tions are weak in the sense that the magnetic field lines 
are only weakly perturbed. The predictions of the LV99 
model including the analytical scaling given by Eq. [S] 
have been successfully testecQ with extensive MHD sim- 
ulations in Kowal et al. (2009). 

It is clear from Eq. ((5]) that the level of turbulence con- 
trols the reconnection rate in the LV99 model. If the tur- 
bulence level is low and the region of magnetic field dif- 
fusion, i.e. field wandering region, is thinner than the A 
given by Ohmic diffusion, the reconnection happens with 
the Sweet-Parker rate. This may explain the accumula- 
tion of magnetic flux prior to Solar flares. However, out- 
flows can destabilize the system by inducing turbulence 
and increasing reconnection/outflow rate (see Lazarian 
& Vishniac 2009). We, however, do not believe that the 
case of marginal initial turbulence is applicable to the 
magnetotail. 

The LV99 model is a model of volume-filled reconnec- 
tion with magnetic filaments/flux tubes filling the wide 

5 The issue of what is happening in 2D reconnection in presence 
of turbulence is rather controversial (see Loreiro et al. 2009, Kulpa- 
Dubel et al. 2010). Unlike 3D reconnection, the X-points must play 
role there if the reconnection is fast (Matthaeus & Lamkin 1986, 
Servidio et al. 2009). However as the nature of turbulence and 
reconnection are different in 2D, we feel that the 2D studies cannot 
clarify much in the physics of the actual astrophysical reconnection. 

6 Weak and strong do not reflect the amplitudes of waves, but the 
interaction between the oppositely moving Alfvenic wave packages. 
For our case the largest amplitude waves are in the regime of weak 
turbulence, but smaller amplitude small scales ripples are in the 
regime of strong turbulence. 

7 Testing of whether reconnection is fast, i.e. independent of 
resistivity, is tricky with present day diffusive codes. Therefore it 
is the successful testing of the analytical predictions that give us 
confidence in the results. 



outflow regioifl The LV99 presented such a model of re- 
connection and observations of the Solar magnetic field 
reconnection support the volume- filled idea (Ciaravella & 
Raimond 2008). In this wide region magnetic field lines 
shrink, converting magnetic energy into other forms. In 
the absence of cosmic rays, the magnetic energy is being 
transfered into kinetic energy of the outflow, but in the 
presence of energetic particles a portion of this energy 
can be channeled into the acceleration of these particles. 

Figure [6] exemplifies the simplest realization of the 
acceleration within the reconnection region expected 
within LV99 model. As a particle bounces back and 
forth between converging magnetic fluxes, it gains en- 
ergy through the first order Fermi acceleration described 
in de Gouveia dal Pino & Lazarian (2003, 2005, hence- 
forth GLOJ0, see also Lazarian 2005). The first order ac- 
celeration of particles entrained on contracting magnetic 
loop can be understood from the Liouville theorem, i.e 
the preservation of the phase volume which includes the 
spatial and momentum coordinates. As in the process of 
reconnection the magnetic tubes are contracting and the 
configuration space presented by magnetic field shrinks, 
the regular increase of the particle's energies is expected. 
The requirement for the process to proceed efficiently is 
to keep the accelerated particles within the contracting 
magnetic loop. This introduces limitations on the parti- 
cle diffusivities perpendicular to magnetic field direction. 
The subtlety of the point above is related to the fact that 
while in the first order Fermi acceleration in shocks mag- 
netic compression is important, the acceleration via LV99 
reconnection process is applicable to incompressible flu- 
ids. Thus, unlike shocks, it is not the entire volume that 
shrinks for the acceleration, but only the volume of the 
magnetic flux tube. Thus high perpendicular diffusion 
of particles may decouple them from the magnetic field. 
Indeed, it is easy to see that while the particles within a 
magnetic flux rope depicted in Figure 6 bounce back and 
forth between the converging mirrors and get acceler- 
ated, if these particles leave the flux rope fast, they may 
start bouncing between the magnetic fields of different 
flux ropes which may sometimes decrease their energy. 
Thus it is important that the particle diffusion parallel 
and perpendicular magnetic field stays different. Particle 
anisotropy which arises from particle preferentially get- 
ting acceleration in terms of the parallel momentum may 
also be important. 

The energy spectrum was derived in GL03 following 
the routine way of dealing with the first order Fermi ac- 
celeration in shocks (see Longair 1992), namely, by con- 
sidering the acceleration rate and the loss rate of ener- 
getic particles without taking into account the backreac- 
tion of the accelerated particles on the flow. This way 
GL03 obtained: 

N(E)dE = C ■ E- 5/2 dE, (6) 
i.e. the spectrum of energetic particles with the spectral 

8 We would like to stress that Figure l5l exemplifies only the first 
moment of reconnection when the fluxes are just brought together. 
As the reconnection develops the volume of thickness A gets filled 
with the reconnected 3D flux ropes moving in the opposite direc- 
tions. 

9 The mechanism published eventually in 2005 paper was made 
public as 2003 arXiv preprint, which motivates us to refer to the 
original on-line publication. 
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index similar to that of galactic cosmic rays (see also its 
derivation in LO09). 

More recently, Drake, Swisdak & Shay (2006, hence- 
forth DSS06) approached the problem of back- reaction of 
accelerated particles on magnetic loops produced by re- 
connectiorPl The authors conjectured that the backre- 
action can be described by the term (1 — 8ire/B 2 ), where 
e is the energy of accelerated particles, and obained the 
spectrum of accelerated particles —3/2 rather than —5/2. 
We feel that the introduction of the backreaction into the 
acceleration process is an important process that requires 
further studies. 

The idea of the acceleration of protons by magnetic re- 
connection has been supported by numerical simulations. 
In Drake et al. (2010) two dimensional calculations of the 
particle acceleration within contracting magnetic loops 
were presented. The first realistic three dimensional cal- 
culations which use the actual reconnection simulations, 
rather than artificially creating loops, were presented in 
Lazarian et al. (2010). More numerical studies of the ac- 
celeration are necessary, as the realistic simulations are 
extremely challenging. 

We note that the acceleration described in GL03 is 
different from electric field acceleration described e.g. in 
Haswell et al. (1992) and Litvinenko (1996). The accel- 
eration discussed there within Sweet-Parker reconnection 
layers is inefficient, first of all, because only a tiny rate of 
magnetic energy release in reconnection processes is dic- 
tated by exceedingly slow reconnection rates. Moreover, 
it is rather problematic to confine the particles within 
the reconnection Sweet-Parker layer. Any field wig gles 
would naturally result in particles leaving the laveira As 
a result the acceleration gets inefficient. 

The estimate of the maximal energy of protons which 
can be accelerated through this process can be obtained 
through the usual arguments that the Larmor radius 
should not be larger than the size of the magnetized re- 
gion L zone (Longair 1997) 

E _, 10l 3 eV (_!_)(_W_), (?) 

which sets the limit of energies, i.e. E max , through ap- 
pealing to the fact that protons of larger energies cannot 
be confined by magnetic fields to experience the acceler- 
ation through multiple bouncing back and forth between 
the reconnecting magnetic fluxeJ^I. 

The properties of the heliotail, especially at large dis- 
tances from the Sun, are not well constrained at the mo- 
ment (see Pogorelov et al. 2009b, Pogorelov 2010, private 

10 Unlike GL03 the study in DSS06 treated particle acceleration 
in 2D reconnection layers, where shrinking filaments degenerate 
into closed loops. The latter study also appeals to the physics of 
collisionless reconnection. However, from the point of view of the 
acceleration process, these differences do not matter much, pro- 
vided that the reconnection processes are fast and the reconnected 
magnetic fields fill the outflow region. The latter condition we feel 
is somewhat more difficult to realize in 2D picture of DSS06, but 
we do not dwell on these details. 

11 It is shown in LV99 that the probability of magnetic field 
return into a reconnection layer is low. Thus the accelerated parti- 
cles are bound to leave the zone of accelerating magnetic field fast 
without getting much of energy gain. 

12 This does not preclude cosmic rays of higher energies to ex- 
perience additional acceleration via electric field as in Litvinenko 
(1996), but we do not discuss this possibility in the paper. 
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Fig. 6. — Cosmic rays spiral about a reconnected magnetic field 
line and bounce back at points A and B. The reconnected regions 
move towards each other with the reconnection velocity Vr. The 
advection of cosmic rays entrained on magnetic field lines happens 
at the outflow velocity, which is in most cases of the order of Va- 
Bouncing at points A and B happens because either of streaming 
instability induced by energetic particles or magnetic turbulence in 
the reconnection region. In reality, the outflow region gets filled in 
by the oppositely moving tubes of reconnected flux which collide 
only to repeat on a smaller scale the pattern of the larger scale 
reconnection. Thus our cartoon also illustrates the particle accel- 
eration taking place at smaller scales. From Lazarian (2005). 

communication), which makes it challenging to identify 
the precise values for magnetic field strength and scales 
involved. A simple estimate of the scale L zone and mag- 
netic field B can be obtain if we accept the solar wind 
velocity 450 km/s (see Pogorelov et al. 2009b) the 11 year 
cycle of magnetic field can create create areas of field of 
one direction of the order of 10 16 cm. Assuming that 
the magnetic field is about 1 /LtG, one can get particles 
with the energy of dozens of TeVs. While the detailed 
calculations, should produce more accurate estimates of 
the magnetic field and the scales involved, we feel that 
we are getting right ballpark numbers. In fact, we can 
predict that, unless some processes of field amplification 
operate in the turbulent heliotail, the acceleration of the 
cosmic rays of energies much larger than 10 Tev is rather 
unlikely by magnetic reconnection. 

6. DISCUSSION 

This paper attempts to explain the cosmic ray excess 
in the range from 50 GeV to 1-10 TeV as arising from 
magnetic reconnection in the magnetotail. The high en- 
ergy cut off observed corresponds roughly to what is ex- 
pected from the reconnection events. Indeed, it is virtu- 
ally impossible to explain the acceleration of higher en- 
ergy particles with the mechanism, unless appealing to 
some hypothetical magnetic field acceleration processes. 
The difference in the distribution of the excess of low and 
higher energy particles in our scenario arises from higher 
efficiencies of scattering for low energy particles. 

Our relation of the observed cosmic ray anisotropies to 
the heliotail is also supported by the sideral daily varia- 
tions of galactic cosmic rays observed with the Tsumeb 
neutron monitor's hourly count during 1977-2000 re- 
ported in Karapetyan (2010). There it was argued that 
the observed cosmic ray excess has heliotail rather than 
galactic origin. No physical explanation of the excess was 
provided there, however. On the contrary, we relate the 
excess of the cosmic rays with the acceleration process 
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induced by magnetic reconnection. 

The paper has an exploratory character, as the quanti- 
tative description of mechanisms of cosmic ray accelera- 
tion in the reconnection regions are at its infancy. Unlike 
shock acceleration, which is the subject of long history 
(Axford et al. 1977, Krymsky 1977, Bell 1978, Blandford 
& Ostriker 1978) and extensive literature (see Gaisser 
1990, Malkov & Diamond 2009). The existing analytical 
models (de Gouveia dal Pino & Lazarian 2003, Drake et 
al. 2006, Drake et al. 2010) are insufficiently elaborated, 
while the numerical simulations (see Drake et al. 2010, 
Lazarian et al. 2010) are rather idealized. Nevertherless, 
the accumulated evidence is suggestive that the process 
of the acceleration can be efficient. 

We note that dealing with the acceleration of protons 
we do not distinguish between the collisionless recon- 
nection and the LV99 model of reconnection. We feel 
that, provided that the reconnection regions are thick 
and filled with the reconnecting shrinking loops, the ac- 
celeration happens the same way for the two cases. We, 
however, claim that the formation of such regions may 
be somewhat problematic within the paradigm of colli- 
sionless reconnection. At the same time, this situation is 
a clear consequence of the LV99 model, which appeals to 
the ubiquitous astrophysical turbulence to enpower it. 

While we accept that the model of acceleration in the 
reconnection regions require more study, we would like to 
stress that the proposed scenario has several attractive 
features. First of all, it allows us to address the issues 
of cosmic ray excess over the entire range of energies 
that these particles were observed. Moreover, the fact 
that the excess of the particles is observed in the direc- 
tion of the magnetotail is suggestive that the processes 
in magnetotail are involved. In addition, the alterna- 
tive mechanisms of producing the excess apparently have 
self-evident problems, as we discussed in the paper. We 
would like to stress, that none of the alternative mech- 
anisms provides a unifying explanation for the existence 
of the excess over a range of energies reported by differ- 
ent groups (see §2). In this situation we think that the 
proposed mechanism should be considered seriously. 

We argue that the localized excesses of cosmic rays in 
the multi-TeV range and the tail-in excess below the TeV 
range are related by the same phenomenology. Within 
our approach the localized regions of the TeV cosmic rays 
are related to the sites of acceleration via reconnection. 
The lower energy particles can be accelerated over ex- 
tended regions of the magnetotail; they are also expected 
to experience more scattering prior to reaching the ob- 
server at the Earth. More elaborate modeling of both 
cosmic ray acceleration in reconnection regions and the 
propagation of cosmic rays in magnetotail should provide 
detailed predictions to be compared with observations. 

Our recent approaches to this problem combine the 
advances in understanding of the statistical structure 
of MHD turbulence, in particular, tensorial structure of 
Alfvenic, slow and fast modes (see Cho, Lazarian & Vish- 
niac 2002, Cho & Lazarian 2002, 2003, Kowal & Lazarian 
2010), analytical description of the propagation of cosmic 
rays (see Yan & Lazarian 2002, 2004, 2008), and test- 
ing of analytical predictions numerically using magnetic 
fields obtained through numerical simulations (sec Bcrcs- 
nyak, Yan & Lazarian 2010). As we described above, 
simulations of particle acceleration in turbulent recon- 



nection is under way. Thus we hope to have in future a 
self-consistent numerically-tested picture of the acceler- 
ation and propagation. Combined with the advances in 
simulations of magnetotail magnetic fields (see Pogorelov 
et al. 2010) this gives hope that we will have a set of 
quantitative detailed predictions for the cosmic ray ex- 
cess arising from the mechanism described in the paper. 

Numerical modeling will also clarify the role of the 
second order Fermi acceleration which arise from turbu- 
lence induced both by magnetic reconnection and exist- 
ing within the magnetotail. The corresponding processes 
have been discussed at length in the literature (see La 
Rosa et al. 2006, Petrosian et al. 2006, Yan et al. 2008), 
but we expect the second order Fermi process to be less 
efficient than the first order Fermi acceleration and there- 
fore to be subdominant. 

In situ measurement of the excess of the energetic par- 
ticle acceleration within reconnection regions could be 
beneficial as well. So far, the attempts of measuring of 
such an excess were not successful (Gosling et al. 2005ab, 
2007, Phan et al. 2007). We believe that this is due to 
the X-point, Petscheck-type reconnection is inefficient in 
the acceleration of cosmic ravJ 1 ^! (see LO09). This type of 
reconnection was sought in the studies above. The mech- 
anism we discussed above (GL03, DX06) appeals to the 
thick extended reconnection regions. Such regions nat- 
urally emerge in the magnetotail as turbulent magnetic 
fields of opposite polarity are being pressed together (see 
Figure 4). 

The GL03 acceleration process has been already em- 
ployed in Lazarian & Opher (2009, henceforth LO09) to 
explain the origin of anomalous cosmic rays, whose mea- 
surements by Voyagers seem to contradict to their origin 
within most of the accepted models of shock accelera- 
tion. A similar conclusion was also obtained in Drake et 
al. (2010) where the process of collisionless, but volume- 
fillin£3 reconnection was discussed. Whether small-scale 
reconnection is collisionless or collisional does not play a 
role for LV99 model. The testing of this fact was suc- 
cessfully performed in Kowal et al. (2009), where the 
plasma effects were simulated through the introduction 
of the anomalous resistivity. 

Unlike LO09, where the acceleration of energetic par- 
ticles in the reversing field of the heliosheath was consid- 
ered, we consider the magnetic field reversals in the he- 
liotail. While the field reversals in the heliosheath arise 
from Sun's rotation with the magnetic axis being not 
parallel to the rotation axis, the reversals in the heliotail 
arise from the 11 year solar cycles. As a result, the scale 
of the reversals is expected to be much larger. This pro- 
vides a possibility of accelerating higher energy particles. 

7. SUMMARY 

In this paper we combine data from different experi- 
ments to prove that there exists a statistically significant 

13 The inefficiency of X-point reconnection arises from both tiny 
amount of magnetic energy being released in the small reconnec- 
tion region and in inefficiency of slow shocks (see Beresnyak, Jones 
& Lazarian (2010) in accelerating energetic particles. An alterna- 
tive explanation of the fact was given by H. Karimabadi (private 
communication) who observed in his kinetic simulations extended 
periods of the stagnation of the X-point reconnection. 

14 We stress this, as the original models of collisionless recon- 
nection (see Shay et al. 1998, 2004), unlike LV99, stressed the 
importance of pointwise, i.e. X-point, reconnection. 
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excess of energetic particles in the direction of the so- 
lar system magnetotail. We proposed an explanation of 
this excess as arising from the acceleration of energetic 
particles in reconnection regions along the magnetotail. 
These regions arise as oppositely directed magnetic field 
of the solar wind is pressed together in the magnetotail. 
The change of the magnetic field polarity arises in our 



scenario due to the well established solar cycles. 
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APPENDIX 

COLLISIONLESS EFFECTS IN MAGNETOTAIL RECONNECTION 

While the LV99 model provides fast reconnection, i.e. the reconnection that does not depend on resistivity, without 
appealing to any collisionless plasma effects, for some of the small-scale events, e.g. for the acceleration of low energy 
electrons the microphysics of reconnection may be important. In the reconnection process described by LV99 model 
the reconnection speed is given by Eq. ([5]) . The same speed can be obtained if one considers the local reconnection of 
flux tubes. It was shown in LV99 that the probability of magnetic field lines which ones entered the reconnection layer 
to reenter the reconnection layer again is low. As a result it was shown in LV99 that the global reconnection rate is 

Vrec, global ~ L / \\\V rec .i oca i (Al) 

where V rec ^ oca i is the velocity of reconnection within the small-scale local Sweet-Parker reconnection events, i.e. 
reconnection events on scale Ay depicted on the lower panel of Figure [5] It is easy to show that assuming that if Aii and 
A_l are related through the GS95 critical balance, namely \\\/Va ~ Aj_/u, Eq. (IA1[) results in too high reconnection 
rates even for Sweet-Parker reconnection at scales A^. The naively obtained rates much exceed those provided by 
Eq. ([5]), which proves that the Ohmic resistivity effects are not the bottleneck of the LV99 model (see more details in 
LV99). 

We will use Eq. (|Alj) ([5]) to establish the scale of the local Sweet-Parker events \\y. 

A,, « LS- 1 ' 3 (L/l) 1/3 (V A /V t ) 4/3 (A2) 
which results in the thickness of the Sweet-Parker layers being 

A turb « LS- 2 ' 3 (L/l) 1/6 {VA/Vif 3 . (A3) 

Those Aturb correspond to X pe rp in the lower panel of Figure O 

The corresponding thickness is much smaller that the thickness of the laminar Sweet-Parker layer, which is LS 1-1 / 2 , 
which makes according to Eq. ([3]) the collisionless effects important for the small scale heliotail reconnection. This, 
however, is not going to change either the rate of magnetic reconnection or the acceleration of protons. The presence 
of collisionless effects can affect the acceleration of electrons, which is, however, not the process that we deal in this 
paper. 

We should add that the issue of fluid being collisionless or collisional is not so simple in the presence of turbulence. 
Collisionless fluids are subject to instabilities, which reduce the effective degree of their collisionality via inducing 
resonant scattering of particles. Such collisions mediated by magnetic field decrease the mean free path of particles. 
For instance, in Lazarian & Beresnyak (2006) the problem of the collisions in a fluid of was treated self-consistently 
in the presence of the gyroresonance instability. It was demonstrated there that turbulent compressions of the fluid 
on the mean free path induce instability at the gyroradius which decrease the mean free path. For thermal plasma 
particles firehose and mirrow instability should also be important (see Schekochihin et al. 2010). As a result, we 
can state that turbulence both decreases the scale over which elementary reconnection events take place, potentially 
allowing the elementary reconnection events to proceed in a collisionless fashion, but at the same the compressions 
arising from the turbulent cascade(3 can decrease the mean free path of the particles. The details of these interesting 
processes are not important for the acceleration of cosmic ray protons that we deal in this paper. 
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